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Abstract 

We derive the loop equations for the d-dimensional n-Hermitian 
matrix model. These are a consequence of the Schwinger-Dyson equa- 
tions of the model. Moreover we show that in leading order of large 
N the loop equations form a closed set. In particular we derive the 
loop equations for the n = 2 k matrix model recently proposed to de- 
scribe the coupling of Two-dimensional quantum gravity to conformal 
matter with c > 1. 
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1 Introduction 



The last three years have seen a rapid development in the knowledge relative 
to two-dimensional quantum gravity coupled to conformal matter with c < 1. 
This has been possible by the use of matrix model techniques and the large 
N expansion, combined with the double scaling limit U Very recently a new 
interpretation of the double scaling limit as a finite size rescaling has emerged 
011 based on the renormalization group applied to large N models |§|J. 

Although the new techniques have provided a very deep understanding 
of two-dimensional quantum gravity coupled to c < 1 conformal matter, a 
barrier has emerged at c = 1, avoiding the continuation of the previous results 
beyond this point. The same barrier appears in the continuum treatment of 
non-critical strings. || 

Various approaches to surpass this problem have been suggested. Among 
them we mention the use of reduced models to discuss the d-dimensional 
matrix model J7| which is the right model to discuss the d-dimensional string, 
Very recently we have obtained the loop equations for this model for arbitrary 
N and show that they form a closed set in leading order of large iV||. For 
the moment they have not been solved but we see them as an important 
exact result that may permit to cross the barrier at c = 1. 

Last year a 2 n -Hermitian matrix model in zero dimensions have been 
proposed to describe a system with c = n/2 coupled to 2d gravity 0. So a 
solution of this model is another way to cross the barrier. 

Four years ago we obtained the loop equations for the zero dimensional 



two- matrix modelflTU], extending the solution of Metha for the partition 



function [11]. Last year we were able to solve these equations completely fl2 



See also 113 



Motivated by these recent progress in the solution of the Two matrix 
model, in the present work we will derive the loop equations for the n- 
Hermitian matrix model in any space-time dimension d and for arbitrary 
N and show that they form a closed set in leading order of large N. They 
are a consequence of the Schwinger- Dyson equations of the model and the 
type of coupling among the different matrices. The couplings we consider 
include the g-state Potts matter fields coupled to Two Dimensional Gravity 
of [0 as a particular case. 

This paper is organized as follows: In section 2 we review the Hidden 
BRST method |14| and use it to derive the loop equations in section 3. In 
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section 4 we discuss in more detail the zero dimensional case. In particular, 
we compare our equations with the two matrix equations of [ 10 1 . In section 
5 we present our conclusions and comment on future work. 



2 Hidden BRST Symmetry and Schwinger- 
Dyson Equations 

The usual method to derive Schwinger- Dyson equations (SDe) involves using 
the invariance of the path integral measure under field translations. Imparl- 
ances of the action itself are not relevant for this derivation. In this section 
we review an alternative derivation which has been proposed some years ago 
Hyp : By making use of a BRST symmetric extension of any action S, we can 
derive all SDe's as BRST supersymmetric Ward identities of the new action. 

To simplify the presentation, consider the path integral describing one 
bosonic field <f(x): 

-s(<p) 



[dip\e- b ^. (1) 

The functional measure [dip] in the last equation is invariant under field 
translations: 

<p(x)^<p(x) + e(x). (2) 
The invariance of the functional measure implies the following statement: 

From (3) we get 

6F __ 5S 

for any F = F(ip). These are the SDe's for the theory. Notice that, in 
general, S is not invariant under (2) and that this lack of invariance plays 
no role in the derivation of the SDe's We now introduce auxiliary fermionic 
variables ip{x), $( x ) an d insert a trivial factor of unity 

- J dxip(x)ip(x) 
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inside the partition function. The resulting action S is invariant under the 
following transformation: 



S[ V ^M = S[<p] + fdx$(x)1>(x) (5) 

5ip(x) = i]){x) (6) 

5ip(x) = (7) 

6$(x) = -5S/5^{x) (8) 



Associated with the BRST-like symmetry, (6)-(8) is a set of Ward identities; 
for example, unbroken BRST implies: 

(5[F(<p)j(y)]) = 0, 

which is 

Computing the average with respect to the fermionic variables we recover eq. 
(4). Notice that the BRST transformation we have just defined commutes 
with any symmetry the original action S may have. In particular, if the 
original action is invariant under a U(N) group, the BRST symmetry will 
commute with this group also and then its implications will be true order by 
order in the 1/N expansion. 

In the next section we use the Schwinger-Dyson BRST symmetry to derive 
the loop equations for the ci-dimensional n-Hermitian matrix model. 

3 Loop Equations 

In this section we will derive the loop equations for the model defined by the 
action: 

S = f d d x\Y, tr{\d kl d^M{ X) k)d»M{x, l)+V l {M{x, £ c tJ trM(x, i)M(x, j)} 

J k,i 1 1 i^j 

(9) 

M(x, I) is an iV x iV Hermitian matrix defined on the (euclidean) space-time 
point x, and Vi is a local function of M(x, I). 
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This action includes the sytems studied in (we follow their notation) if 
we take: 

d TT i = — — — ff ' exp Bt ■ t' (10) 

(2sinh2/3)" FH v ' 

T . r' = aia[ + a 2 cr 2 . . . + a n a' n (11) 

r = a 1 a 2 ...cr n (12) 

c TT i = d TT i (13) 

where a = 1, 2, . . . , q, for the q-state Potts matter fields coupled to gravity. 
The extended action and the BRST transformation are: 

S ext =S + Jd d xEitrtlj(xJ)tlj(xJ) (14) 
8M(x,l) =^{x,l) (15) 

6ip(x,l) =0 (16) 

Introduce the following representation of M(x,l): 

N 

M{x, l) = J2 m «( x > l ) Ta { x , 0> ( 18 ) 

a=l 

where the projectors of the matrix M(x, I) satisfy: 

JV 

Q = l 

trT a (x,l)T p (x,l) =5 af3 (20) 
T a {x,l)TP{x,l) =5 a pT a (xJ). (21) 

If 5 denotes the BRST variation (see last section), we get the following BRST 
transformation for the projectors T a and eigenvalues m a (see the Appendix 
of I): 

§T a ( X I) =T T a (x,l)^x,l)T' } (x,l)+Tl } (x,lMx,l)T a (x,l) 

8m a {x,l) =ti^{x,l)T a {x,l). (23) 



4 



Consider the following functional of X(x, I) 



N 



Kij = [ d d x £ £ A Q (x, l)T a (x, (24) 

^ J Q=l 

The basic object we will use to write the loop equations is: 

«[A] = e lK (25) 

It fullfils the identity: 

J U l X \, = f 1 dte ltK tT a (x, l)e^ K , (26) 
dX a (x, I) Jo 

which implies: 

6tTu[X] - ttrT a (xJ)u[X] (27) 



5X a (x,l) 

The trace is taken with respect to the internal indices only. From u[X] we 
can compute all symmetric combinations of products of T a . For instance: 

5 2 u[M | TfeQr 2 , f 2 )T^(m, h) + Tft (xx, /i)T fe (x 2 , f 2 ) 

U=0 — ~ 



SX f3l (x l ,l 1 )5Xf 32 (x 2 j2) 

(28) 

The Schwinger- Dyson equation is: 

(5tr[u[X]$(x,l)\) = (29) 

i.e. 

(*r[<J«[A]^(x, Z) + *r[u[A]<ty(x, /)]) = 0. (30) 

Since 

6rj>(x,l) = Y,d lk DM(x,k) + Y,c lj M(x,j) - A t (M) (31) 

where 

A l (M) = V l '(M(x,l)). (32) 

We get: 

truKdj*CIM(a;,A;) - A^M) + ^cijM(x,j)} = 

k m 

Y^[^2(dik(m a (x, k)U x + 2d fl m a (x, k)d„ + n x m a (x, k)) - A a (x, k)5i k )truT a (x, k) + 

a k 

cijm a (x, j)truT a (x, j)} .(33) 



We have that: 

truT^l) = -i^^. (34) 

and 

N 



5u 



Jq I a=l (3^a 



m a (x, I) — mp(x, I) 
We compute the fermionic average to get: 

(5u$Cx,l)) = (iV f 1 Xa [ Xl l ), - A/3( f ; l \ tru[t\]T p (x, l)tru[(l-t)\]T a (x, I)). 

(36) 

Since T a (x,l) appears only once in the trace, we can express it as a 
derivative of u: 

fra[(1 _ t)A]n , ^^fcM (37) 

We get the following loop equation: 

ri ^ A a (x, Z) - A /3 (x, /) 1 d>tru[£X] 5tru[(l — i)A] . _ 
m a (x,l) -mp(x } l)t{l-t) 5\p(x, I) 5\ a (x,l) 

- (E E ( x > fi jrr^ + 

Y^ d ik(m a (x, k)n x + 2d^m a (x, k)d^ + a x m a (x, k) - A a (x, k)6i k ) f^^j.s ])>( 38 ) 

which is valid to any order in the 1/N expansion. 

We can simplify the equation if we restrict ourselves to the leading order 
in 1/N. Then we can apply Witten's factorization property |15|t the loop 
equations to prove that: 

F(M) ~ (F(M)) + 0(N~ a ), a > 
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for any £7(iV)-invariant F. In particular, we obtain that: 

m a (x,l) ~ (m a (x,l)) = (m a (0,l)). 

That is, in leading order of large N the loop equation becomes a closed 
system for u[X]: 



X a (x, I) - \ p (x, I) 1 8tru[tX] 5tru[(l - t)X] 

(JjL~ 



^ a Jo m a (l)-mp(l) t(l - 1) 5Xp(x,l) 5X a (x,l) 

5tru[X] 

E<k(m.,(fc)Q, - AJk)S tl ) ^% ], (39) 

We could derive this closed system because we chose to write the Schwinger- 
Dyson (loop) equations in terms of T a (x, I) rather than M(x, I). Also by us- 
ing the BRST Schwinger- Dyson transformation we could avoid this difficult 
change of variables in the path integral and the (highly) non-trivial measure 
in the (m a , T a ) basis. 

The loop equations must be solved subject to the following boundary 
conditions: 

tru\ x =o = N 
Stru 



5X a (x,l) 

5 2 tru 



5X a (x,l)Xp(x,l) 



o = -5«p (40) 



u = constant (A independent) is a particular solution of the loop equa- 
tions, but it does not satisfy the boundary conditions. We must look for a 
non-trivial solution. 

We get an equivalent formulation of the loop equations noticing that only 
the derivative of tru appears in the loop equations. So identify: 

< 41 > 
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So that the loop equation is equivalent to the following system of equations: 

f 1 dt K(x l) -X,(x l) 1 [tx] { l)v[{1 _ t)xu z) = 
m a (l)-mp(l) t(l-t) 

£ di k (m a (k)D x - A a (k)5 lk )v[X] a (x, k)] (42) 



<Ju[A] a (x,Z) _ (^[A]^?/, j) 



(43) 



8Xp(yJ) 5\ a (x,l) 

v[\] a (x,l)\^ = i (44) 

fa[A] a (g,Q 
xw n U=o = -<W (45) 

Notice that: 

trM(x 1 ,h)M(x 2 ,l 2 )...M(x n ,l n ) = ]T m ai ...m an (l n )trT ai (x 1 ,h)...T an (x n ,l n ) 

a\...a n 

= — m ai ...m an [trT ai (a;i, Zi)..T a "(x n , /„) + all permutations of indices a] (46) 

ai...a n 

Then it is enough to know tru[\] to get all correlations of the matrix 
M(x, I) that are symmetric under permutations of the space-time points 

{Xii l{) ■ 

A continuum version of the large-iV loop equation is obtained by the 
usual manipulations. That is we introduce the functions X(z, x, I) and m(z, I) 
defined for < z < 1 and the density of eigenvalues pi, and make the 
following identifications: 

\ a (x,l)=\(^,x,l) (47) 
a 



m«,(() = v / JVm(-,() (48) 
Mm,) = £ (49) 
J dmipi(mi) = I. (50) 
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The large- N loop equation becomes: 

P [ dm 1 dm 2 pi(m 1 )pi(m 2 ) [ dt ^— v -^Il K m 2,y,l) 1 x 

J Jo mi — m-i t{± — t) 

6tru[tX] 6tru[(l-t)X] ^ f , , N tru\X] 

TV7 77 — r— 7 — = — > Cu / dmpAm)m-— r 

5X(m 2 ,y,l) 5X{m x ,y,l) j£ h J HA 1 SX(m,y,j) 



dmpi(m)(ma - Ai(m))— -; (51) 

J oX{m,y,l) 



P stands for the principal value of the integral over m,. 

4 The Zero dimensional case 

In this section we discuss the zero dimensional case. The equations are: 



E fdt 



X a {l) - X P {1) 1 Stru[tX] 6tru[{l - t)X] 



o m a (l) - m p {l) t(l - t) 5X/3(l) 5X a (l) 

-aE<.M)||-A.wS]. (52) 



We want to check that we reproduce the equations of JTIJ for n = 2. In this 
case we have: 

u[X] = 1 + i £ X a (l) + ^[£(A Q (1) 2 + A Q (2) 2 + 

a,l 1 a 

2 A a (l)A /9 (2)tr(T a (l)r a (2))] A a (l) 3 + A a (2) 3 + 

3£(A Q (1)%(2) + A Q (l)A / ,(2) 2 )tr(T Q (l)T /3 (2)) + o(A 4 ) (53) 

Inserting this into the loop equations we readily get the system of UlOH . In the 
same way we can get the equations for the higher correlators of T a (l), Tg(2). 

It is, of course, a very interesting problem to see if the methods of |L2| 
can be extended to the n-matrix chain. 
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5 Discussion and Open Problems 



We have been able to derive the loop equations for the d-dimensional n-matrix 
model. They are a consequence of the Schwinger-Dyson equations satisfied 
by the Green functions of the model. Due to the factorization property of 
U(N)-invariant operators in large N the loop equations form a closed set that 
may be the starting point for a non-perturbative description of the system. 
In particular our result applies to the q-state Potts matter fields coupled to 
two dimensional gravity that have been previously proposed to describe the 
non-critical string with c > 1. 

There are several interesting open problems. It is known that the loop 



equations for d — models can be realized as Virasoro|l| and W 3 JT 
constraints on the free energy of the system. We wonder whether this is true 
for the present case too. In addition the loop equations could be used as 
the starting point for a different approximation to the Physics of the model, 
perhaps on the lines suggested by the solution of the Two matrix model in 

0- 
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